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ABSTRACT 

", In fully general relativity, we calculate the images of the radiatively inefficient accretion flow (RIAF) 

■ surrounding a Kerr black hole with arbitrary spins, inclination angles, and observational wavelengths. 

' For the same initial conditions, such as the fixed accretion rate, it is found that the intrinsic size and 

, radiation intensity of the images become larger, but the images become more compact in the inner 

region, while the size of the black hole shadow decreases with the increase of the black hole spin. With 
the increase of the inclination angles, the shapes of the black hole shadows change and become smaller, 
even disappear at all due to the obscuration by the thick disks. For median inclination angles, the 
radial velocity observed at infinity is larger because of both the rotation and radial motion of the fluid 
f***- ' in the disk, which results in the luminous part of the images is much brighter. For larger inclination 

angles, such as the disk is edge on, the emission becomes dimmer at longer observational wavelengths 
(such as at 7.0mm and 3.5mm wavelengths), or brighter at shorter observational wavelengths (such 
as at 1.3 mm wavelength) than that of the face on case, except for the high spin and high inclination 
images. These complex behaviors are due to the combination of the Lorentz boosting effect and the 
radiative absorption in the disk. We hope our results are helpful to determine the spin parameter of the 
black hole in low luminosity sources, such as the Galactic center. A primary analysis by comparison 
with the observed sizes of Sgr A* at millimeters strongly suggests that the disk around the central 
black hole at Sgr A* is highly inclined or the central black hole is rotating fast. 
Subject headings: black hole physics - Galaxy: center - radiation transfer - sub-millimeter 
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More and more observations indicate that supermassive black holes reside in most local g alaxies, and the co- 
evolution of galaxies, quasars and supermassiv e black holes was realized in last decade (e.g., iFerrarese fc Merrittl 
120001 ; [Gebhardt et alj|2000t iHopkins et all 2008). From the stellar dynamics, it is not surprising to find that the com- 
| pact radio source Sagittarius (Sgr) A* at the center of ou r Galaxy is associated with a~4x 10 6 M Q massive dark object 
. (jSchodel et al.l 12003: Gh ez et all 120051 iGhez et all 120081). A natural candidate of the massive dark object is massive 
black hole, as emphasized bv lKormendv fc Richstond (|1995). they argued convincingly that this object is a black hole, 
and the other alternatives to a black hole should be ruled out. Some of these alternati ves have alre ady been excluded. 
They include a cluster of dark objects, such as the ste llar remnants or brown dwarfs (!Maoz|[l998f). a binary made of 
point masses, self gravitationa lly supported Fermi ball (Miller 2006; [Tsiklauri fc Viollierll 19981 : iMunvaneza et aT1ll998l : 
lYuan. Naravan. fe Rccs 20 04)) ■ Another alternative which is still wait ing for being excluded is self gravitationally 
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supported boson ball ([Torres et al.l 12000: Yu an. Naravan. fe Re cs 2004). However, the problem of this model is that 



the mass of bosons is arbitrarily ass umed to ma ke a massive boson ball, and the energy dissipation between bosons 
may lead to the collapse of the ball (Maoz 1998). 

Recent very-long baseline interferometry (VLBI) observation at the wavelength of 3.5mm shows its intrinsic size is 
0.126±0.017 mas, about ~ 1AU, and the resulting lower limit of its mass density is 6.5 x 10 21 M Q pc~ 3 (|Shen et al.l 
120051 ). and the VLBI observation at 1.3mm gives the intrinsic size 0.037+q'q^q mas, and the resulting lower limit of 
its mass density is 9.3 x 10 22 M & pc~ 3 (jDoeleman et al. 2008), which strongly supports that Sgr A* is a supermassive 
black hole. The milli meter VLBI observation has been used to estimate the parameters of SgrA* accretion flow 
(jBroderick et al.ll2008h . In the near future, higher r esolution observations of Sgr A* at infrared (iPaumard et aT1l2005h 
and sub-millimeter wavelengths will be available (jDoeleman fc Bower! 2004; M ivoshi et ail 12004 iFish et al.l I2008L 



which can provide a direct test of the physical processes under the strong gravity, such as the dynamics of the accreted 
gas, complex trajectories of photons, a nd so on. Therefore, it is necessary to investigate the images of the flux and 
polarization of realistic accretion flows ([Yuan et al.ll2004l iBroderick fc Loebll2006al : iHuang et alJl2007L l2008f), and the 
images and l ight curves associated with a hot spot in the accretion flow near the black hole horizon (jBroderick fc Loe"bl 
120051 l2006bf ). the hot spot is proposed to explain the observations of near-infrared and X-ray flaring of Sgr A*. 

It is generally believed that Sgr A* was a low luminosity AGN whose activity was switched off due to insufficient 
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gas supply (|Naravanll2002t ). Multi- wavelength observations put strict constraints on theoretical accretion models of 
the massive black hole at the Galactic center. Among the theoretical models, the advection-dominated accretion flow 
(ADAF) model, also called the radiatively inefficient accretion f low (RIAF) model is claimed to be very su c cessfu l 
in modeling the multi-wavelength spectrum from radio to X-ray (jNaravan et al.lfl995l 119981 : lYuan et alJl2003L 12004). 
Based on the RIAF model, the images of Sgr A* at millimeter have been calculated, the size of the images is c ompared 
with observations to put an independent constraint on the RIAF model (jYuan et al.1120061 iHuang et al.ll2007t ). and the 
shap e of the images, especially the location of the image centroid of Sgr A* can be used to determine the black hole 
spin (|Broderick fc Loebll2006ah . 

There are several shortcomings in the previous calculations: first, although the treatment of the photon trajectories 
is within the framework of fully general relativity by using a ray tracing method, the accretion model is Newtonian, 
which means the global structure of the accretion flow is based on the Newtonian dynamics. The general relativity 
effect for the accretion flow is mimicked with Paczynski & Wiita potential. Therefore, the radial velocity of the 
accretion flow very close to the black hole is larger than the speed of light and thus is unphysical. Second, the massive 
black hole is assumed to be a non-rotating (Schwarzschild) hole in I Yuan et al.l (|2004h . 

In this paper, following lManmo"tol (|2000f ). the global structure of the relativistic ADAF around a Kerr black hole is 
derived to calculate the images of Sgr A* surrounding a Kerr black hole with an arbitrary black hole spin and viewing 
angle at several millimeter wavebands. Our treatment is within the framework of fully general relativity. A brief 
introduction to relativistic ADAF model is given in §2. The ray tracing method and the radiation transfer in curved 
space-time are discussed in §3. In §4, we summarize the key steps for calculating the images of ADAF of Sgr A*. The 
results and discussions are presented in §5 and §6, respectively. 

2. GENERAL RELATIVISTIC RIAF MODEL 
In the Boyer-Lindquist coordinates (t, r, 9, </>), the Kerr metric is 

ds 2 = -e 2v dt 2 + e 2 ^{d<f) - ujdt) 2 + e 2 ^dr 2 + e 2 ^d0 2 , (1) 

where 

2 ,_SA 2 „_ Asin 2 9 S 2ar 

e A' E , A' ' A' 1 ' 

A = r 2 -2r + a 2 , S = r 2 + a 2 cos 2 9, A = (r 2 + a 2 ) 2 - a 2 A sin 2 9. (3) 

Here, the notation in the seminal paper (jBardeen et al.lfl972f ) has been taken, and the geometrical units G — c — 1 
are used throughout this paper. The accretion flow has both azimuthal and radial motion in the accretion flow. The 
velocity field of the fluid is described by fl and V, where fl is the azimuthal angular velocity of the fluid and V is 
the physical velocity of the fluid wit h respect to the co-rotat i ng reference frame which is carried by the observers with 
the fo ur velocity u» = u*(l,0,0,fi) ([Abramowicz et al.lll997l : IPeitz fc Applll997t iGammie & Pophamlll998t IManmotol 
2000). The four velocity of the fluid is (see the Appendix for the details) 

u M = {lr^e~\ j r r e-^ , 0, u*fi). (4) 

In order to solve the global solution of the RIAF in Kerr metric, we follow IManmotol ([2000) to formulate the fully 
relativistic hydrodynamical equations for RIAF. For simplicity, they assume the height of the accretion flow H < r, 
therefore the Boyer-Lindquist coordinates (t,r,9,(f)) are expanded around the equatorial plane up to (z/r) terms to 
obtain a cylindrical coordinates (t, r,(j>, z). Another approximation is generally taken in previous investigations: the 
vertical integration. T herefore, all the ph ysical quantities are taken their values at the equatorial plane. We employ the 
approach suggested by Manmotc] (|2000f ) to calculate the global structure of an accretion flow surrounding a massive 
black hole in the general relativistic frame, which allows us to calculate the structure of an accretion flow surrounding 
either a spinning or a non-spi nning black hole . All the radiation processes are included in the calculations of the global 
accretion flow structure (see IManmotol 120001 for the details and the references therein). Integrating these equations 
from the outer boundary of the flow at r = r ou t inwards the black hole, we can obtain the global structure of the 
accretion flow passing the sonic point smoothly to the black hole horizon. Here, we allow the accretion rate of the 
accretion flow to vary with radius due to the putative winds (see [Blandford & Bcgc lmanlll999L for the details), i.e., 
we assume the ac cretion rate M(r) = Mq{t /r on t) s , where Mq is the accretion rate at the outer radius of the ADAF 
(lYuan et al.M2003h . 

The global solution of general relativistic RIAF (GR-RIAF hereafter ) has been obtained in the previous works, and 
the e f fects of black hol e rotation have been discussed in the great details (|Gammie fc Pophamll9 98: Poph am fc Gammid 
1998; Manmoto 2000). In order to help understand our numerical results, the global solution of ADAF is shown in 
FigHJ The basic parameters are taken as follows: the black hole mass of Sgr A* = 4 x 10 6 M Q , the wind index 
s = 0.25, the fraction of the dissipated power directly heating the electrons 5 = 0.3, the accretion rate at the outer 
radius Mq = 0.7 x 10 _6 A/©yr~ 1 , the ratio of the gas pr essure to the total pressure (3 — 0.95. These parameters are 
the same as the RIAF parameters in I Yuan et al.l ([20061 ) 's work, except the slightly difference in accretion rate (see 
SJS] for the discussion). It should be emphasized that the ion and electron temperatures increase significantly within 
10 2 i? g with the increase of the black hole spin parameter from a = to a = 0.998, and decrease, with less difference, 
with its decrease from a — to a = —0.998. Similar changes appear in electron number density and magnetic field 
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Fig. 1. — Global solution of GR-ADAF. The solid lines represent the results for a = 0, the dotted lines the results for a = 0.998, the 
dotted-dotted-dashed lines the results for a = 0.5, the dashed lines the results for a = —0.5, and the dotted-dashed lines the results for 
a = —0.998. Except the upper left panel, the lines above the solid ones represent the results for prograde disks (a > 0), while the lines 
below the results for retrograde disks (a < 0). Upper left: Radial velocity V as a function of radii. The lines below the solid ones are for 
a > 0, while the lines above are for a < 0. Upper middle: Coordinate angular velocity Q. Upper right: The half height of the disk H. 
Lower left: The number density of ions versus radii. Lower middle: The temperature of ions and electrons. Lower right: The magnetic 
field B which is deduced from the global solution of ADAF. 

strength. Because the radio emission is mainly from the synchrotron radiation which depends sensitively on electron 
temperature, the effects of the black hole spin on the radio emission is mainly through its effects on the electron 
temperature. 

3. RELATIVISTIC RADIATION TRANSFER OF PHOTONS FROM GR-RIAF 

In this paper, the images of Sgr A* modeled as a GR-RIAF surrounding a massive black hole are calculated. In our 
calculations, we do not consider the Compton effect of photons at radio wavelength, because the relativistic electrons 
will boost the low energy photons to higher energy which is beyond our interest. As a bunch of photons travel along 
their geodesies, some of them will be absorbed due to synchrotron self-absorption in the accretion flow, while some of 
new photons will be added due to the emission of the flow. To solve the radiation transfer of radio photons, we should 
first determine their trajectories using the ray tracing method. 

3.1. Ray Tracing Method 

In the Kerr space-time, all of the four constants for the motion of a test particle have been found by iCarterl (|1968l) . 
they are the energy at infinity 

E = - Pt , 



the effective angular momentum, 

the Carter constant, 

and the Hamiltonian constant, 



Q 



a 2 E 2 cos" 



: P>p, 



Li cot" 0, 



H=\e= 



(5) 
(6) 
(7) 

(8) 
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where e = for massless particles such as photons, or e = 1 for massive particles. Therefore, given three constants 
of motion, i.e. E,L Z ,Q, the trajectory of a photon can be determined. As the trajectory of a photon is independent 
on its energy, it can be completely de t ermin ed by the two dimensionless parameters A = L z /E and q = Q 1 / 2 /E. 
Furthermore, Cunningham & Bardecn ( 1973) pointed out that the two parameters A and Q are related to the two 
impact parameters a and [3 that describe the apparent position of the image on the sky as seen by an observer at 
infinity who receives that light ray, 

VpW\ A 



a = - TT = -- -IT- > ( 9 ) 

V P ( ) ) r^oo Sin ^obs 

= (^ij") = + « 2 COs2 d obs ~ A 2 cot 2 obs ) 1 / 2 = Pohs , (10) 

where o bs is the angular coordinate of the observer at infinity and is the four momentum in the locally non-rotating 
frame (see the Appendix). The direction of (3 points along the projected direction of the spin axis of black hole. Due 
to the azimuthal symmetry, the photon trajectory originating from a position (r, 0) to the observer at infinity whose 
angular coordinate is o b s can be solved by the following integral equation: 

f°° dr f e ° hs dQ 

± / = ± / = P, (11) 

Jr ^R{r;\ iq ) Je ^0(0; A, q) 

where 

R(r; A, q) = r 4 + (a 2 - A 2 - q)r 2 + 2[q + (A - af]r - a 2 q, (12) 

6(6»; A, q) = (q + a 2 cos 2 - A 2 cot 2 0). (13) 

The signs in the above integrals must be the same as those of dr and d9 to guarantee the integral P is always 
positive and increasing along a photon's trajectory. Therefore, the signs change at the turning points in r and 0, 
where R(r;X,q) — or 0(0; A, g) = respectively. It is important to realize that both R(r) and 0(0) are quartic 
polynomials of r and respectively, therefore both integrals in Eq. (|11[) can be integrated in terms of elliptic functions 
(jChandrasekharlll983l ; iRauch fc Blandfordl[l99l iCadez et al.1ll998f ). 

Owing to above observations, a photon trajectory can be determined by using a ray tracing method: first, given 
a and (3 which determine the position of the image in the photo-graphical plate of an observer with o bs, the two 
constants of motion are obtained, according to Eqs. (|9"l fTU|) . Second, given the integral P, which can be understood as 
the affine parameter along the photon's trajectory, the position of the photon can be determined using the expression 
of Eq. fTTj) as elliptic functions: 

r = r(P;X,q), (14) 
= 0(P;A,g). (15) 

3.2. Relativistic Radiation Transfer 
Given the two constants of motion of a photon, its momentum has components, 

dx v ( V~R i 

= 9nv-j— = (kt,k r ,k e ,k<),) = I -1,±— ,±VQ, A j (16) 

where a is the affine parameter of photon's trajectory, and E is absorbed into a. The signs in the above equation are 
the same as those of dr and dO. Suppose the photon frequency observed by the observer at infinity is fobs, then the 
emitting frequency at the locally rest frame of the fluid (LRF) is 

fern = -(k^U^Vobs, (17) 

where is the four velocity of the fluid. In the LRF, the radiation transfer equation is (Rybicki & Light manlll979f) 

dl(v em ) 



dL 



-K(v em )I(i> em ) +r)(v em ), (18) 



where / is the radiation intensity of photons and dL is the element of physical distance of photon traveling in LRF. 
As T = I /v 3 is a Lorentz invariant, the above equation can be re-written as follows, 

— — — = - K (jy em )J(i/ cm ) -) — . (19) 

aL v em 

Since the speed of light is c — 1, we have 

dL dx^ dx^ en-, . , . , 

— = — — = — ef LRF = -k^, (20) 
da da da p 
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Therefore, the radiation transfer equation in curved space-time reads ([Fuerst fc Wd 12004?) . 

■ = -kfj,u^{-K{v em )l{v eia ) + r](i' em )/i' em ] (21) 



da 

= [-V e mK(v eIa )Z{v<im) + ^cmV^cml/^obs- (22) 

It is well known that T,vk{v), and r](i/)/i/ 2 are Lorentz invariants (jRvbicki fc Lightmanl 119791 ) . therefore, the above 
radiation transfer equation is written in the Lorentz invariant form as we expected. 
Since 

dv i — 

S— =±vR, (23) 
du 

therefore, 

da = -EdP (25) 
Introducing k = u em K(u em ) / v bs and rj = ^(t / e m)/(j / em l 'obs)) an d the formal solution to the transfer equation reads, 

J(oo) = J expf-y Tt(P' )EdP" j rj{P )EdP' , (26) 

where P max is the maximal trajectory integral tracing from the infinity to the innermost emitting position. 

4. IMAGES OF GR-RIAF AT MILLIMETER WAVELENGTHS 

After solving the global structure of the RIAF around a Schwarzschild or Kerr black hole, the physical quantities at 
the equatorial plane have been obtained, such as: 

n e (f, z = 0), m(r, z = 0), V r (f, z = 0), fi(f, z = 0). (27) 

In the following calculations, the isothermal approximation is taken in the vertical direction, therefore, the density 
distribution in the vertical direction is: 

n e (r, z) — n e (f, z = 0) exp(— z 2 /2H 2 ), 
m(r,z)=m{f,z = 0)exp(-z 2 /2P 2 ), 

V r (f,z) = V r (f,z = 0), 
Q(f,z) = n(f,z = 0). 

At any position along the photon trajectory, such as (r(P),9(P)), we can determine m m (P) and fc M (P), so we get 
^em(P) = {—u^{P)k tJj {P))v \ is . The magnitude of magnetic filed B(f, z) is determined by the local physical conditions: 

P gas (f, z)=m(r, z)k B Ti + u e (f, z)k B T e , 

Pmag(r,z)=Pgas(r,z)(l-/3)//3, 
B(r, Z) = V^Pmag 

The emissivity and absorption coefficients at each position can be calculated, considering the bremsstrahlung and 
synchrotron radiation mechanism. The Compton scattering effect is ignored, because it generally produces photons at 
higher wavelength band which are not interested in this work. Finally, the radiation intensity received by the observer 
at infinity can be calculated according to Eq. (|26|) . 



and 



5. NUMERICAL RESULTS 

In the work of lYuan et~aTI (|2006f ). the size of Sgr A* was calculated in the RIAF model of lYuan et~aTl (|2003h . In this 
RIAF model, the accretion rate is assumed to be M{r) = Mo(r/ r out ) s due to the presence of winds, where Mo is the 
accretion rate at the outer radius of the RIAF (jYuan et al.ll2003l ). The black hole mass o f Sgr A * M b h = 4x 10 6 M Q is 
adopted in their work, which is different from Afbh = 2.5 x 10 6 M Q adopted in lYuan eF al. (2003J). Thus, they tuned the 
parameters of the RIAF in order to reproduce the observed spectrum in different wavebands, which leads to s = 0.25, 
the fraction of the dissipated power directly heating the electrons 8 — 0.3. Their calculations of the RIAF structure 
were carried out by adopting the Pseudo-Newtonian potential. Compared with the results in fully general relativistic 
frame in the case of a Schwarzschild black hole (i.e., a = 0), their global solution of the RIAF deviates only slightly 
in region out of ~ 20P 5 , but to a certain extent i n the inner region within ~ 20R g . In this work, we simply adopt 
the same RIAF parameters as in lYuan et all ([2006) in all our calculations independent of black hole spin parameter a, 
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Fig. 2. — Intensity of the ADAF around Kerr black holes with different spins at the wavelengths of A = 7.0mm, 3.5mm, 1.3mm (columns 
from left to right) as a function of a (or (3) when O ^ S = 0°. Results for different spin parameter a are shown in line styles the same as 
those in Fig[T] The horizontal lines show the intrinsic size of the disk which is defined according to Eq. (28)- (29). 
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Fig. 3. — Images at wavelength of 7.0mm of the ADAF of Kerr black holes with spins a = —0.998, —0.5.0, 0.5, 0.998 (from left to right) 
and with the viewing ang les 6> obs = 90°,45°,0 (from top to bottom). The spin of the black hole is along the perpendicular direction. The 
diameter of the dash-dotted circle is twice the photon capture impact parameter (2\/27M). The diameter of the dashed circle is the intrinsic 
size (-D c h) of each image, see definition in the text. To show relative radiation intensity, the range of the intensity is scaled i together for 
all the panels. 

except that a slightly lower accretion rate (M = 0.7 x 10 -6 M Q yr -1 , 70 percent of that adopted by I Yuan et alj (|2006l )) 
at the outer boundary of the RIAF is adopted in our calcula t ions, which can compensate the discrepancy between 
our GR-RIAF solution and the RIAF solution of lYuan et al.l (|2006| ) in the inner region of the flow near the black 
hole with a — 0. We find that our global GR-RIAF solution can reproduce the broad spectrum from radio to X-ray 
quite well for a non-rotating black hole while viewed face-on. We treat this as our fiducial model. In principle, we 
should have to tune the RIAF parameters to fit the observed spectrum, the size of the scattered images, and even the 
po larization properties, w henever the black hole spin parameter a and the inclination angle o bs are specified as done 
in IBroderick et aTl (|2008D . but this is very cumbersome and beyond the scope of this work, which should be carried 
out in our later work. To focus on the influence of black hole spin on the general properties of the images at millimeter 
wavelengths, as a first step, we will not change any RIAF parameters in our calculations for the GR-RIAF structure, 
except for different values of black hole spin parameters. Because the non-thermal electrons are artificially added to 
explain the observed long wavelength (> 7mm) radio spectrum, we consider the radio emission is contributed only 
from the thermal electrons in the following calculations. 

Following the procedures described in the previous sections, we calculate the images of GR-RIAF around a black 
hole with arbitrary spins, inclination and at different observational wavelengths. Our results are shown in Figf2][Hl 
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Fig. 4.— As in FigE] but images at the wavelength of 3.5mm. 
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Fig. 5. — As in Fig|3] but images at the wavelength of 1.3mm. 



respectively. 

Figure [2] shows the results for face on case, i.e., the radiation intensity of GR-RIAF around a black hole as a 
function of a (or /?, because of the symmetry of the system) for a = 0, ±0.5, ±0.998 and for millimeter wavelengths, 
A bs = 7.0, 3.5, 1.3 mm, in which we are most interested. Here we define the intrinsic size -D c h of the disk at a given 



TABLE 1 

Size (diameter) of the shadow of the black hole 



Spin of the hole 


A = 7.0mm 


A = 3.5mm 


A = 1.3mm 


(a) 


(M) 


(M) 


(M/fias) 


0.998 


2.8 


2.8 


2.8/13.8 


0.5 


5.2 


5.2 


5.2/25.6 





5.6 


5.6 


5.6/27.6 


-0.5 


(5.0 


6.0 


6.4/31.5 


-0.998 


6.4 


6.8 


7.6/37.5 



frequency as 

f I u {a,f3)dn ohs =e / ' I v (a,f3)dQ ohs , (28) 

Jp<D ch /2 J 

where p = ^/ a 2 + (3 2 is the distance to the mass center of the hole, £l bs is the observed solid angle, and the coefficient 
is the ratio of emission within the full width of half maximum (FWHM) of a 2-dimensional(2-d) Gaussian function 
G(a,0) to its total emission, i.e., 

_ Jp<FWHM/2 ^( a ' 0)d&obB 

£ ~ ~ /G(a,/3)^ obs ' 

By convention of VLBI observation, FWHM is used to measure the size of a 2-d Gaussian distributed source. 

For an image with the standard Gaussian distribution, the intrinsic size D c h defined here reduces to its FWHM, 
therefore, for an image deviated from the Gaussian distribution, the intrinsic size can be considered as its equivalent 
FWHM. As seen from Fig[2l the radiation intensity I u increases with the increase of BH spin ranging from a — —0.998 
to a = 0.998, while the intrinsic size of the disk, marked as region with radii covered by the horizontal lines for 
each case, increases in the same way. These results are expected from the effect of BH spin on the global solution 
of GR-RIAF (see FigUJ): the larger the BH spin, the higher the electron/ion number density and their temperature, 
and the larger magnetic field in the disk. Because the cmissivity at radio wavebands is dominated by synchrotron 
radiation, larger spin leads to higher radiation intensity of the disk. As has been found in previous investigations, 
roughly the global solutions of GR-RIAF for a > di ffer from that for a = more than the corresponding a < cases 
(|Gammie fc Po pham 1998; Pop ham fe Gammid fl998) . This tendency is also reflected in Figj2]for longer wavelength 
A bs = 7.0,3.5mm. With the decrease of the wavelength, the emitting region becomes smaller, because the physical 
conditions in the inner region change dramatically, consequently, the difference among the final intensity in the inner 
region becomes larger, especially for a > cases. The effect of the observational wavelength on the intrinsic size of 
the images is also clearly shown in FigO the shorter the observational wavelength, the smaller the emitting region of 
the disk. 

For a Schwarzschild black hole, the maximal impact factor for photon capture is V27M, therefore, the typical size of 
the BH shadow is the circle with diameter of 2V27M. The size of the black hole shadow, region without any emission, 
are listed in Table [TJ In the last column with A = 1.3mm, the corresponding angular size is also listed. 

In our calculation, the accretion flow extends to the horizon of the hole . It is well known that the larger the 
spin of the hole, the smaller its horizon and marginally stable orbit r ms . This tendency is clearly shown in Table [TJ it 
is found that the shadows of black hole decreases with the increase of spin of the hole, if the disk is prograde with the 
hole, whilst, if the disk is retrograde with the hole, the size of the shadow does not change significantly. For the same 
spin of the hole, the siz e of the shadow at different wavelength does not change significantly (see also IFalcke et alJ 
2000; Brodcr ick fe Loebll2006al |b1: lTakahashill2004f ). 

Figure [3]l5] shows the images of the GR-RIAF around a black hole with different spin a, the inclination angle of the ob- 
server at infinity # bs, and at the different wavelength A b s - The observed images at 7.0mm for the prograde/retrograde 
disks with the black hole are shown in Fig[31 Each column corresponds to the results for different spin of the hole, 
from left to right, the spin of the hole is a — —0.998,-0.5,0,0.5,0.998, respectively. Each rows corresponds to the 
different 6* b s , from top to bottom, 6* b s = 90°,45°,0°. For each image, its intrinsic size is plotted in a white dashed 
circle, emission within which is equal to e which is defined in Eq. (29) times its total emission. The photon capture 
region for a Schwarzschild black hole is also plotted in a white dash-dotted circle for reference. All images in this figure 
are scaled in the same way to show clearly the relative radiation intensity for different spin and inclination. 

In Figure [31 for the same inclination, the effects of spin on the size and intensity of the disk are similar: the larger 
the BH spin, the larger the emitting region, i.e., the intrinsic size, and the higher the radiation intensity. As increasing 
the inclination angle, the shadow of the images decreases due to the obscuration of the GR-RIAF with a finite height 
which is about H < r. For the extreme case of 8 bs = 90°, there is no shadow at all. When d ^ s = 45°, the luminous 
part of the images which is located at the left and top (for a > cases) , or right and top (for a < cases) of the image 
is brighter than the luminous part of the images for o b s = 0°. This is due to the Lorentz boosting effects: at the left 
and top (or right and top) of the disk, the fluid has the larger radial velocity to the observer at infinity because of the 
rotation and radial motion of the accretion flow. For o bs = 90°, the radio emission is separated into two arcs for high 
spin (e.g. a > 0.9), while the weak radio emission appears between two arcs for low spin. Our analysis is as follows: 
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the photons from the central part travel through the most dense region near the equatorial plane of the disk, which 
results in the great absorption of the emission, therefore, the resulting emission is much weak. Radiative absorption is 
the reason why the luminous parts of the images are not at the central regions near the equatorial plane of the disk, 
and the intensity of the luminous parts is also lower than that of the corresponding face on cases. As discussed above, 
the intrinsic size decreases with the increasing of the inclination angle. However, such decrease becomes less significant 
for larger spin. For the a — 0.998 case, the intrinsic size nearly does not change. This can be explained as follows: for 
larger inclination angle, the obscuration of the disk decreases the emission, as well as the size of the emitting region. 
On the other hand, the light bends more significantly, bringing emission from the more inner region, therefore both the 
emission and emitting region increase. For an extremely-rotating (a = 0.998) black hole, the emitting region increased 
by its extremely curved space-time can almost compensate that decreased by inclination, although the emission is not 
compensated enough. 

The images at 3.5mm and 1.3mm are shown in Fig. |4] and Fig.[5j respectively. The behavior of the images at 3.5mm 
and 1.3mm with varying the spin and the inclination are very similar to those at 7.0mm. The main differences are 
reflected in the intrinsic sizes of images. At shorter wavelength, the disk appears to be smaller than those at longer 
wavelength, this is because the shorter wavelength photons are from the more inner parts of the disk. For shorter 
wavelength photons emitting from the inner parts of the disk, their optical depth might be so significant, causing the 
shadow in the images of the disk even when 9 Q h s = 90° for low spin. The size of the shadow is similar to the size of 
the circle with a radius of the photon capture impact parameter (\/27)M. When the disk is edge on, for the longer 
wavelength of the observed photons, because the optical depth of the emitted photon near the equatorial plane is 
not so large that the shadow does not appear. For A b s = 1.3mm and median spin a < 0.5, the luminous parts of 
the images are brighter than these of the corresponding face on cases, so that the larger inclination results in higher 
emission. This is because the Lorentz boosting effect dominates the radiative absorption. Exceptionally, at high spin, 
such as a = 0.998 and # b s = 90°, the emitting region is near to the horizon of the black hole, thus, the radiative 
absorption dominates the Lorentz boosting, as a result, the intensity of the luminous parts is lower than that of the 
corresponding face on cases as in Fig. [3] and Fig. 0] for longer wavelengths. 

6. DISCUSSIONS 

In this work, the images of the GR-RIAF around a Kerr black hole are simulated. The main results are summarized 
as follows: 1) the effects of BH spin: the larger the spin, the smaller the si ze of black hole s hadow, which is similar 
to the previous results (e.g., iFalcke et al.1 120001 iBroderick fc Loebl 2006a, b; Takahashi 2004)). The larger the spin, 
the larger the intrinsic size of the image, and the brighter the whole disk, especially its inner part. These effects are 
more signific ant in a prograde disk th an a retrograde disk. We should emphasize that t his contrasts with the result 
obtained in (|Broderick fc Loebl l2006aH bT). The reason is that IBroderick fc Loebl ( 2006a, b) assumed a parameterized 
ADAF solution for a spinning black hole which is required to fit the observed spectrum, therefore, they found that the 
thermal electron density must decrease with increasing spin, as well with increasing inclination. More important is 
that the non-thermal electron density must increase, in order to produce the radio spectrum. However, in this work, 
we do all the calculation in the light of a fiducial GR-ADAF model in which the accretion rate is fixed. Consequently, 
the electron density increase with increasing spin (see Fig.[T]). 2) the effects of the inclination: the larger the inclination 
angle s, the smaller the size of black hole shadow, and the smaller the intrinsic size of the image. This is consistent 
with (Noble et al. 2007). If the disk is edge on, the shadow are even obscured at all. The central part of the image is 
both darkened by the absorption and brightened by the light bending, so that how the brightness of the image changes 
depends on the observational wavelength. 3) the effects of the observational wavelength: the shorter the observational 
wavelength, the smaller the intrinsic size of the image. 

Although a larger spin makes the intrinsic size enlarged, it makes the total emission enhanced to a greater extent. 
This tendency differs from face-on to edge-on cases. At 7.0mm, for instance, as shown in Fig{3j an extremely-rotating 
black hole (a = 0.998) enhances > 60% total emission than a non-rotating black hole (a = 0) does. For the face-on 
case, however, the intrinsic size is enlarged only by ~ 15%. For the same comparison at 3.5mm, as shown in FigSl 
the total emission is enhanced by w 99% while the intrinsic size is enlarged by only ss 23%. Therefore, more emission 
comes from the inner pa rt of the image if t he black hole spin increases. 

As discussed in §5, in I Yuan et al.l (|2006l ), the global structure of the RIAF were calculated by adopting Pseudo- 
Newtonian potential, and taking the non-thermal electrons into consideration. The disk is viewed as facing on, and 
the thickness of the disk is ignored . It was found th at the scattered images are more extended than the observational 
images at both 7.0mm and 3.5mm (jShen et al.H2005h . In our GR-RIAF model with the same parameters, the scattered 
images are even bigger. This is because the lower electron temperature and density are obtained in GR-RIAF model, 
enlarging the emission region. However, the overestimated intrinsic sizes do not imply an anti-rotating black hole, 
which predicts smaller intrinsic sizes. As mentioned above, this is because the anti-rotating black hole also predicts 
much lower total emission. The required emitting region should be enlarged a lot to produce the observed emission. 
On the other hand, a co-rotating black hole is preferred, since the required emitting region could be more compact. 

For a primary application, we adjust the index of the wind index s and energy fraction of non-thermal electrons 
r\ to reproduce the spectrum of quiescent state of Sgr A*. Only face-on cas e is investigated a s an example. For the 
Schwarzschild case, we choose s = 0.285 and rj — 0.4%, different from those in I Yuan et alj (|2006f ) due to the relativistic 
effects. In order to produce the observed the spectrum both in radio and sub-millimeter band, s should be decreased, 
while i] should be increased for a faster co-rotating black hole. Contrarily, for a faster anti-rotating black hole, s 
should be increased while rj should be decreased for the same purpose. In FigJBl the predicted spectra, the FWHM 



10 



36 



^ 35 : 



a 0.99.S 
a =0.5 
a=0 
a =-0.5 
a =-0.998 



X 7mm 



—i 



34 : 



60 

5 33 



32 






0.80 




0.78 


(maj 


0.76 


S 




WH 


74 


to 






0.72 




0.70 




0.55 




0.50 


"a 
'i 






0.45 


X 








to 






0.40 




0.35 



10 



11 12 13 
log 10 (v/Hz) 

A 1.3mm 



14 



15 



X 3.5mm 



d 

u 

T3 



T3 



o 

o 




JCMT - CARMA 

JCMT - ARO/SMT 



0.35 
„ 0.30 

to 

0.20 

0.15 
0.30 

_ 0.25 

a 

'g 

2 0.20 
to 

0.15 
0.10 



Fig. 6 

Sgr A*, 



1 2 3 4 -l.o -0.5 o 0.5 l.o 

Baseline (G 1) a 

. — Spectra, full-width of half maximum of the scatter-broadened image at 7.0mm and 3.5mm, and visibility profiles at 1.3mm of 
predicted by face-on GR-RIAF with a = -0.998, 0.5, 0, 0.5 and 0.998. See the details in the text. 



of the images at 7.0mm and 3.5mm, and the visibility profiles at 1.3mm are shown for five examples with different 
spins, with the interstellar scattering effect having been taken into consideration. As expected, the smaller sizes are 
obtained at all th e three wavelengths for a higher spin. Compared with the measured sizes at 7.0mm and 3.5mm by 
iShen et all (|2005f ). the predicted FWHMs are plotted in the solid and dashed lines for averages and limits in the right 
two panels of Fig[6j respectively. It is strongly indicated in Figj6]that a fast co-rotating black hole is preferred, if the 
disk is facing on. Notice that the FWHM of the minor axis at 7.0mm overestimates the observation, even with an 
extremely rotating black hole. This may imply an inclined disk w hich makes the intrinsic sizes smaller (see section 5). 
By analyzing the existing spectral and polariztion measurements, iBroderick et al.l (|2008f ) found that low values of the 
inclination angle are excluded, which is consistent with our analysis. 

The VLBI observations also predict that the black hole shadow structure is hoped to be resolved if the observational 
wavelength decreases to 1.3mm. The scattering size extrapolated from previous longer wavelength observations is 
~ 23/ias, which means any structure with size larger than it would not be washed out. As shown in the bottom-left 
panel of FigO the long baseline is predicted for the existence of the null point, which would be an evidence of the 
shadow region of a black hole with high spin. The hi gh spin causes the emi ssion region more compact. Compared to 
the latest visibility detections at 1.3mm reported bv lDoeleman et ail |2008), only two of the five cases, with a = 0.5 
and 0.998, marginally reproduce the correlated density at ~ 3.5GA detected by the JCMT — ARO/SMT baseline. 
These results also suggest that the black hole is fast rotating or the inclination angle is large, which is consistent with 
what suggested by the previous observations at longer wavelengths. As shown in table [1] a fast rotating black hole with 
a > 0.5 predicts size of shadow less than the scattering size, so that its shadow structure is hard to resolved, especially 
for a highly-inclined disk. If the Galactic central black hole is really fast rotating, the future shorter wavelengths 
(sub-millimeter) arc necessary for the detections. 
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APPENDIX 



THE REFERENCE FRAMES 

In the ADAF model, the accretion flow has both azimuthal and radial motion. Following Peitz & Appl (1997) (see 
also Gammie & Popham 1998, Manmoto etal 2000), we introduce three reference frames to show how to describe the 
velocity fields of accreted gas. The first is the locally non-rotating frame (LNRF) which is also called zero angular 
momentum observer (ZAMO), an orthonormal tetrad basis carried by observers whose four velocity in Boyer-Lindquist 
coordinates is = e~ u (l, 0, 0, ui). The four contra-variant basis vectors of LNRF compose a matrix as follows, 



(a) 



(LNRF) : 



e~ v 
e^ 1 



(Al) 



while the four covariant basis vectors compose a matrix as follows, 

e v 

s (LNRF) 





e^ 2 

we* e* 



(A2) 



As well known, the four rows of the matrix give four vectors in the curved space-time, whilst the four columns give 
four vectors in the local inertia frame, i.e. four vectors in the 'Minkovski' space-time. 

The second frame is the so-call co-rotating frame, another orthonormal tetrad basis carried by observers co-rotates 
with the fluid with the coordinate angular velocity f2 — u^/u 1 , i.e. whose four velocity is = 0,0,O). The 
relative physical azimuthal velocity of CRF with respect to LNRF is — e^ _ly (f2 — uj), therefore, the four basis 
vectors of CRF can be obtained by an azimuthal Lorentz boost from these of LNRF, 
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The final frame is the local rest frame (LRF) of the fluid, another orthonormal tetrad basis carried by observers 
moving with the fluid, i.e. whose four velocity is w M = (u*, u r , 0, ifi). Suppose the relative physical radial velocity of 



LRF with respect to CRF is 
boost from these of CRF, 



V, therefore, the four basis vectors of LRF can be obtained by a radial Lorentz 
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